Theory of Raman scattering from Leggett's collective mode in a multiband 

superconductor: Application to MgB2 
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In 1966 Leggett used a two-band superconductor to show that a new collective mode could exist 
at low temperatures, corresponding to a counter-flow of the superconducting condensates in each 
band. Here, the theory of electronic Raman scattering in a superconductor by Klein and Dierker 
(1984) is extended to a multiband superconductor. Raman scattering creates particle/hole pairs. In 
the relevant Aig symmetry, the attraction that produces pairing necessarily couples excitations of 
superconducting pairs to these p/h excitations. In the Appendix it is shown that for zero wave vector 
transfer q this coupling modifies the Raman response and makes the long-range Coulomb correction 
null. The 2-band result is applied to MgB2 where this coupling activates Leggett's collective mode. 
His simple limiting case is obtained when the interband attractive potential is decreased to a value 
well below that given by LDA theory. The peak from Leggett's mode is studied as the potential 
is increased through the theoretical value: With realistic MgB2 parameters, the peak broadens 
through decay into the continuum above the smaller (tt band) superconducting gap. Finite q effects 
are also taken into account, yielding a Raman peak that agrees well in energy with the experimental 
result by Blumberg et el. (2007) . This approach is also applied to the ? = 0, 2-band model of the 
Fe-pnictides considered by Chubukov et al. (2009). 

PACS numbers: 74.70.Ad, 74.70.-b, 74.25.Gz 



I. INTRODUCTION 

In a superconductor, a longitudinal density oscillation 
corresponds to an oscillation of the phase of the super- 
conducting order parameter. As shown by Bogolyubov 
et el} and Anderson,^- this oscillation would amount to 
a collective mode with zero energy at zero wavevector in 
the absence of the long-range Coulomb interaction be- 
tween the conduction electrons. The presence of this in- 
teraction, however, causes the condensate density oscil- 
lation to become the electronic plasma oscillation, typi- 
cally with an energy of 5-10 eV. Leggett showed that for 
a simple model of a two-band superconductor it could 
be possible for the superconducting condensates in each 
band to oscillate longitudinally relative to one another, 
equivalent to an oscillation in the relative phase of the 
two condensates, with the net displacement of the total 
condensate being zerO)^ The frequency of this mode and 
its dispersion with wavevector q would depend substan- 
tially on the effective interband pairing potential. If this 
potential were sufficiently small, but not zero, the small 
q mode frequency would be less than the smaller of the 
superconducting gaps 2A, for the two bands. In the low 
temperature limit this mode would become a true collec- 
tive mode with an infinite lifetime. 

There is good evidence that MgB2 is a multiband su- 
perconductor. Whereas theory suggests that there are 
2 bands primarily of 2pa character and 2 bands of 2p7r 
character,^ most experimental superconducting proper- 
ties can be described as though there were one a band 
and one tt band^^ The a band is nearly two-dimensional. 
Its barrel-shaped Fermi surface is almost cylindrically- 
symmetric about the c axis and the c axis component 



of the Fermi velocity is very small. On the other hand, 
the TT band Fermi surface is three dimensional. On it the 
root mean square Fermi velocity components are large 
and nearly isotropic. 

Electronic Raman scattering measurements in E2g 
symmetry {xx — yy, or xy) show two features due to su- 
perconductivity: a relatively sharp peak at 2Ao- = 13.5 
meV associated with the a band and an onset without a 
peak at 2A7r = 4.6 meV associated with the tt bandfi As 
shown by Blumberg et el., the situation is different for 
Raman in the fully-symmetric Aig symmetry (xx + yy).-^ 
Apart from two-phonon features, an electronic Raman 
peak was found at 9.4 meV, midway between the onset 
of the TT band continuum at 2An^ and the a band peak 
at 2Acr. Because it can decay into the tt band contin- 
uum, the mode signified by this peak cannot be a true 
collective mode. We will show that it represents a res- 
onance mode that follows by continuity from Leggett's 
collective mode when allowance is made for (1) a vertex 
correction that is related to the mechanism that allows 
electronic Raman scattering to couple to it, (2) damping 
due to decay into the tt band continuum, (3) necessary 
kinematic corrections that depend upon q • Vp, where q 
is the transferred wavevector and Wi? is the Fermi ve- 
locity vector, (4) integration over values of q that are 
ill-defined due to optical absorption, all using a realistic 
set of physical parameters that apply to MgB2. 

The present work is based on a paper by Klein and 
Dierker (K-D)^ that uses a Green's function approach to 
describe the electronic Raman response function as re- 
sulting from Kawabata's 4-vertex function^'' The latter 
may be collapsed into a 2-vertex polarization 'bubble.' 
The bare vertex 7^ comes from the electron-photon inter- 
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action and is the sum of two terms: the term, where 
A is the vector potential, couples to the electron den- 
sity operator in first order, and the p • A term couples 
in second order, (p is the electron momentum operator) 
to a density-like quantity through a virtual intermediate 
state. The sum of these terms, 7^, gives the amplitude 
for excitation from near the FS to near the FS of an elec- 
tron from wavevector k to wavevector k + q while the 
incident laser (L) photon is transformed into the scat- 
tered (S) photon with respective wavevector and energy 
transfers q = ki^ — kg and w — uji^ — ujg. For small val- 
ues of q relative to the Fermi wavevector, the vertex 7^ 
does not depend on q. Due to the p • A term, will 
depend on lol when the excitation energy to the inter- 
mediate state comes into resonance with lul- Such res- 
onance Raman effects are probably not relevant for the 
Blumberg et el. data taken with 1.65 eV laser light: Op- 
tical data^ and band structure calculations^ show that 
interband optically-allowed transitions to and from the 
Fermi surface in MgB2 occur at energies well above 1.65 
eV. The 'mass approximation' for the relevant Aig sym- 
metry may then be applied. It says that the electronic 
Raman vertex 7^ is Tq, the classical radius of the elec- 
tron, multiplied by (fixxik) + A*yy(^))j where is the 
inverse effective mass tensor at wavevector k near the FS. 
The holelike a band and the electronlike tt band possess 
^J'xx + ^J■yy values that have opposite signs. 

To allow for the possibility of electronic resonance(s) 
with intermediate state(s), we assume in our derivation 
that 7fc is complex due to the appearance of an imaginary 
term in the energy denominator(s) proportional to the 
inverse lifetime of the intermediate state(s). 



needs. We assume that k, hereafter denoted by fc, is a 
wave vector near the Fermi surface, referring to one or 
more bands, and consider a fc-dependent superconducting 
gap Afe , normal state energy minus the Fermi energy Sk , 
and Fermi velocity Vk — ^ k^k Define 
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The correct branches of the square root and arcsin 
functions are obtained by adding an infinitesimal posi- 
tive imaginary part to the frequency variable oj, which 
is itself assumed to be positive, and by assuming that 
branch cuts are along the negative real axis of the argu- 
ments of the two functions. Note that when q — 0^ pk 
and rfc equal fk, but when q ^ 0^ they no longer equal 

fk- 

Assume zero temperature and that a non-retarded ef- 
fective potential I?^^],/ (pairing correction) acts between 

a p/p pair at k and one at k' and that a potential -D^"]., 
(vertex correction) acts between a p/h pair at k and one 
at k' . We obtain the two vertex equations, rearranged 
from K-D Eqs. (10a, b), that connect the 'dressed' ver- 

tices and for the p/h and p/p-phase channels, 
respectively with the bare p/h vertex 7^: 



II. MODEL AND ITS EQUATIONS 

To carry out the calculation, we will set up the problem 
in a more general fashion and then make the assumption 
that 7fc is a positive constant for k in the tt band and 
a negative constant for k in the a band. These bare 
vertices create a positive charge-density (particle/hole or 
p/h) excitation of wave- vector q in the tt band and a neg- 
ative charge-density excitation of wave- vector q in the a 
band. In the Nambu formalismii used by K-D, these 
p/h excitations are associated with the Nambu ma- 
trix. A generalized vertex equation necessarily couples 
the (p/h) Ta channel with the particle/particle (p/p) T2 
(phase) channel that describes excitations of the conden- 
sate. This is the first of four final state corrections to 
be consider in this paper, the 'pairing correction.' The 
opposite signs of the bare 7^ in the p/h channel produce 
opposite forcing terms in the p/p-phase channel for the 
2 bands and thus enable coupling to Leggett's collective 
mode. The dynamics of the mode are altered because of 
the second final state correction that occurs in the p/h 
channel, the 'vertex correction.' 

We begin by quoting key results from Sections III. 2 
and III. 3 of K-D, somewhat rewritten to fit our present 
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Here dSk = {2'K)^^d^k5{ek) and the Dirac delta func- 
tion 5{£k) forces the k integral to be over the Fermi sur- 
face. Lk = sinh^^(ajfe/Afc) is a log-like term involving 
a cut-off energy ujk and the gap. It also occurs in the 
BCS-like equation for the gap, K-D Eq. (18): 



A, 
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Also needed is the equation for the polarization bub- 
ble, or Raman response function, rearranged from K-D 
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Eq.(lOc): 



B. 



2AI 



(2) 



(7) 

Here the notation B^' ^r-yi^^, q) nieans that the bare 
vertices are jk and 7^, with 7^ dressed by vertex cor- 
rections. The argument {oj,q) refers the fact that fk, 

(2 3) 

Pk, Tfc, and rj^. ' depend on the same argument, omitted 
here for simpUcity. The factor of 2 outside the integral 
is from a summation over spin inherent in the Nambu 
formahsm. 

In general, vertex corrections result from an infinite 
sum of ladder diagrams. This has the consequence that 
the result is symmetric in the two vertices. Thua^^ 



BY^r^{uj,q) = B^j^/(a;,q). 



(8) 



Next we make the third final state correction by cor- 
recting the response function for the long-range Coulomb 
interaction, K-D Eq. (17), to obtain the screened re- 
sponse function 



B^.,A7(w, q) = B^.,r^(a;, q) + 



1 • , 47re^ 
and with Vq = — — 
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Here Vc is the volume of a unit cell. For values of q that 
are small compared to the inverse lattice constant 
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The appearance of 1 on the right side of these equa- 
tions signifies that the relevant bare vertex is unity, 1 , for 
all values of k. On the right side of Ea.(ITO| the second 
term, the so-called 'screening correction,' projects off, in 
a q and uj dependent fashion, from the bare vertex 7^ 
that portion parallel to 1. When 7^ oc 1, the projection 
cancels the first term, and the screened response function 



III. DEALING WITH THE SPREAD IN WAVE 
VECTOR DUE TO OPTICAL ABSORPTION 

The photon wavevectors are proportional to the com- 
plex index of refraction. Thus, for Raman scattering from 
metals, the wavevector transfer is a complex quantity 
q = q' + iq" ■ The experiments of Blumberg et alfi- were 
such that q was parallel to the hexagonal symmetry axis. 
If z denotes the distance along the axis measured from 
the sample surface, then the perturbation that drives 
the Raman response has spatial dependence e'*^' 
This has a Fourier coefficient (i/[-\/27r(g' -|- iq" — q)]. 
The optical penetration depth, {2q")~^, is the order of 



1/lOth the wavelength of light and thus much larger than 
a lattice constant. Therefore, Raman is essentially a 
bulk probe, and excitations with different wavevectors 
q are independent, that is, non-interfering. The physi- 
cal screened response function is obtained by integrating 
i3A(7),7-(w,(j)|l/[\/27r(q' -|- iq" - gf)]p over q from -00 
to 00. The Raman intensity /(w) is minus the imaginary 
part of the result. Using the fact that i?7».A {to, q) is even 
in q, we obtain a result equivalent to K-D's Eq. (4a): 



Im[g^>,A^(c^ + zO+, q)]{q'^ + q"^ + q^)dq 
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IV. APPLICATION TO THE MULTI-BAND 
CASE: SIMPLIFYING ASSUMPTIONS 

We introduce j as a band index and kj to denote 
wavevector values in band j. The superconducting gap 
Aj and the magnitude of the Fermi velocity Vj are as- 
sumed to be constants in each band. Let u denote the 
cosine of the angle between the direction of q and the 
direction of the Fermi velocity v^j at kj. The quanti- 
ties appearing in Eqs.([T]|31) wiU depend on the discrete 
variable j and the continuous variable u: 



Pk I3j{uj,q,u) 
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:(13) 



Pk Pj{uj,q,u) = fj{uj,q,u) 



rk rj{uj,q,u) 



up — (qvju)'^ 



(14) 



(15) 



For n bands, the effective pairing potential is assumed 
to be represented by a symmetric n x n matrix J)^^^ with 
elements Z?^^l when k is in band 7 and k' is in band 7'. 

Let the symmetric n x n matrix D^"-* denote the effective 
potential in the p/h channel responsible for vertex cor- 
rections and allow for the possibility that D*-^-' ^ D^^\ 
With these assumptions, the integrals in Eqs.(j4][7l) sim- 
plify into sums over band index j of integrals over the 
Fermi surface of that band, where the only variable in 
the integrand is u. The integral equations (|4]) and ([5]) 
turn into discrete equations for the dressed vertices 



.(2,3) 



.(2,3) 



(16) 



Upon introduction of the density of states pj for band 
j, the integrals appearing on the right side of Eqs.(|3]- 
[7]) reduce to pj multiplied by an average over u. Let 
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Fj{uj,q), Pj{uj,q), and Qj{uj,q) denote these averages of 
q, u), Pj{uj, q, u), and rj{uj, q, u), respectively. Make 
n vectors F{u},q), P{u},q), and Q{uj,q) out of them, as 
well as n vectors out of pj, Aj, Lj^ bare vertices 7^ and 



dressed vertices Tj ' '(uj,q). For simplicity, we will often 
drop the explicit dependence on {LO,q). Let / denote the 
71 dimensional unit matrix. The dot • denotes the scalar 
product of two vectors or matrix multiplication. 

The existence of Leggett's mode requires the existence 
of the inversei^ 



(17) 



Above, dependence on {i-^,q) remains understood, 
and IIQpII, ||fp|| and ||Pp|| denote n x n diagonal 
matrices with components Qj{uj,q)pj, etc., on the 
diagonal. The symmetric matrix 



U = -A {WA - (VK • A) /) , 



or, in component form 



(25) 



The following matrix equations were obtained. The 
explanation immediately follows. 
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Equations ([TS)) and ([2T|) are matrix versions of Eqs.Q 
and (O. Eg nations and (1^ were obtained by left 
multiplying the matrix equations we obtain from Eqs.([S]) 
and (ini) by Above, when we write, e.g., LpT^'^\ we 
refer to the n-vector with components LjPjT^ . 

Using the result ip = -A^^W^ • A from Eq. ^ 
to eliminate Lp from the first term on the right side of 
Eq.([T^. we solve for F'-'^' in terms of F^'^-'. After inserting 
the result in Eg. dTSl) and solving for r*-^-' (oj, g) in terms 
of 7, we obtain from Ea.(pT|): 



By,r^(c.,g)=7'.il(f^^')(c.,g).7. 



(22) 



Here R^^'^\uj,q) is the response-function matrix in the 
presence of vertex corrections due to pairing (p) and the 
vertex correction (w) due to the residual interaction in 
the particle/hole channel. 



imprints on the matrices R^-^\uj,q) and R^^'^^ (oj,q) the 
effect of the vertex correction due to pairing. U has the 
property that the elements in a given row or column sum 
to zero. Recalling 1, the n-vector consisting of I's. this 
property implies U ■ 1 — 1 ■ U = 0. 

Our fourth correction for final state interactions will 
result by summing bubble diagrams consisting of the re- 
sponse function matrix R^-^'^\uj, q) connected by a sym- 
metric interaction matrix that we will denote by D^''^ ^ 
For clarity, we consider the possibility that D'-^-' 7^ 
^ The series R^P^''\uj,q) + R^P-''\u,q) ■ 

may be summed. The result, to- 
gether with Ea. (|23p . gives a new response function matrix 

fi^P'^'^H^, 9) obeying: 



R 



(p,i',6) 



(27) 



The screened response function Xii'^il) = 
B-y'^A-yi^jq) then becomes 

with R = R^P^''-''\LU,q). 

As discussed in the Appendix, the screening correction, 
the second term, vanishes as q — 0. 



R^P'''\oJ,q) = 



R^PHi^,q) 



(23) 



THE 2 BAND CASE WITH ZERO WAVE 
VECTOR 



where ii'-^-' (w, q), is the response function matrix with a 
correction only due to pairing: 



\Qp\\ + \\Fp\\ 



rU-\\Pp\\ 



\Fp\\ 



(24) 



The remaining discussion concerns the case of n=2. 
Using the above assumption that D^^^ has an inverse 



(and thus the determinant D^^^ 



-{o'f^)^ is nonzero). 



we write the q = Q result for 2 bands as follows. (It can 
also be written in the form of the expression given in 
Eqs. 3, 5 and 6 of Blumberg et al.f, but generalized to 
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include the bubble diagrams.) 
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(29) 
(30) 

(31) 

(32) 

(33) 
(34) 



The screening correction gives zero, as it must for q = 
0. See the Appendix. 

It is a convenient, but incorrect, assumption that final 
state interactions can be neglected. Almost all authors do 
not make any corrections: pairing, vertex, or bubble, rep- 



•){p) 

these quantities are set equal to zero, we obtain 



resented above by nonzero -D^j^ D\j' , and Z?^^^^ — When 



X°{L^,q = 0) = 



-2|72-7iP 

Pi Jl + P2 J 2 



(35) 



This is identical to the result that would be obtained by 
using the expression with screening, Eq.psp. but with 
r(P'^M(^^^ 0) replaced by the bare response function ma- 
trix R°{u},0) = -2 ||F(w,0)p||. The difference between 
X-f{^,q = 0) and x'^i^iQ — 0) is both quantitative and 
qualitative. 



VI. APPLICATION TO FE-PNICTIDES 

Whereas the q — assumption is not realistic for 
MgB2, it is for the Fe-pnictides, magnetic superconduc- 
tors with a layered structure that makes them nearly two- 
dimensional, with very small components of the Fermi 
velocities perpendicular to the layers. For them, we ap- 
ply Eqs. ipSj) through ([M)) to the two-band model with 
extended s-wave symmetry considered by Chubukov et 
ali^. Those authors assume that A2 = — Ai — —A and 
Pi = P2 = P- In their notation, our matrices describing 
pairing, vertex, and bubble corrections would read: 



Ui U3 

Ui Ul 
Ul Ui 



We then find 
X^{cJ,q = 0) = — 

1 - I Ueff 
Ueff = 2ui — U2 — Ui. 



Ui U2 
U2 Ui 



-|72-7i|W(5i) 



(36) 



(37) 



Using renormalization group arguments, Chubukov et 
al. asserted that Ueff is dominated by ui and is positive. 
This means that among all final state interactions, the 
dominant one is repulsive, and it couples p/h pairs in 
one band to those in the other band. They assigned the 



value 0.4 to u 



eff- 



If 



is neglected with respect 



to Ueff, our result becomes similar to theirs: 



X^{uj,q^ 0) 



172 - 7i 



'Pf{jk 



1 - "e/// [JS, 



(38) 



When plotted, Imx7(a; + «0+, 0) for Ueff = 0.4 is similar 
to the plot shown in their Fig. lb. There is a difference. 
Based on an argument involving an energy-dependent 
cut-off of the effective interaction, they used / (^) — 1 in 
place of our / (^). This replacement would give incor- 
rect results if applied to the Ideal Limit to be discussed 
next. 



VII. TWO-BAND CASE WITH FINITE WAVE 
VECTOR IN THE IDEAL LIMIT 

Following Leggett^ we consider the case where the ma- 
terial parameters and regions of interest are such that uj 
and qvj are both assumed to be non-zero but small com- 
pared with Aj. This leads to the approximations 



Fj — Qj 



1 ; P. 



where 



1 



FS 



(39) 



(40) 



The last equality results from the assumption that the 
Fermi surface is spherical. With these assumptions, 
Eqs.(l27]and[24l) lead to the result 



R 



-2 



(41) 

Recall that ||p ^|| and ||c^p|| denote diagonal matrices. 
This gives the following results from Eq.(| 



\^2-li\^H + vW)/Vo 
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(42c) 
(42d) 
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[Recall Eas. ((33l) and (l34l) ]. We derive from the denom- 
inator of Eq. (|42ap the existence of a collective mode at 
an energy lu — u)L{q) obeying 



2 2 



(43) 



with an energy ('mass') oJq and a dispersion velocity Vq. 
These are generalizations of the results given by Leggett^ 
in his Eqs.(3.39). If our parameter Vo were given by 
Vo = + p^^, and taking into account the fact that 
our density of states pi is for one spin, whereas his is 
for both spins, we would obtain his result. Thus, setting 
D*-^' = 0, i.e., the absence of a vertex correction 
in the particle/hole channel and of a bubble correction, 
gives his results for the mode frequency and its dispersion 
with q. 

Use of Eqs. fTT]) and (|42ap gives us an analytical result 
for the Raman intensity in this ideal case. If < Wq, 

/ideoi(w) =0. If W > LOo-, 



lideali^^) ~ 

<i>(g',g",c.)^ 



(44a) 



5/,,16g"c^^$(g',g",a;) 

(Pl + P2)VoUo\/^ 

^/2 _^ y-2g^2 

(44b) 
|72 -7iP(Pi +P2) 



4q" 



Name 


Dii 


D22 


D12 


Pi 


P2 




Value 


-0.47 


-0.1 


-0.08 


2.04 


2.78 




Units 


Ry-c 


Ry-c 


Ry-c 


/Ky-spm-c 


/Hy- Spin- c 




Name 


Ai 


A2 


Vl 




q' 


tt 

q 


Value 


6.75 


2.3 


0.38 


4.7 


1.73 


4.18 


Units 


mcV 


meV 


lO^m/s 


lOV/s 


loV-i 


ioV-1 



(44c) 



TABLE I: Values of input parameters used in the calculation 
for MgB2. Here 'Ry' denotes the Rydberg, and 'c' means 
'unit cell.' To convert the units of the product of q' with v\ 
into meV, multiply by ft = 0.658 ■ 10~^^ meV s. 



On the other hand, for pairing and for the vertex correc- 
tion, no such symmetry arguments apply, and we simply 
assume that D^") = D^^^ = D. 

The material parameters used by Blumberg ei al^ are 
those shown inTable I. The values of g" and are esti- 
mates based on optical data. 4 The other parameters are 
from Liu et al.^. Using these parameters we calculate 
cjo = 12.42 mcV and Vo 0.23 -lO^m/s. This value 
of Wo is much too large for the conditions of the ideal 
limit to be satisfied. Note that roughly Wo oc \J—D\2 
and that Vo depends very weakly on —D12 as long as 
— 2Di2 ^ + P2^ — Dii — D22, which is the case here. 
For example, if —D12 were to equal 0.001 Ry-cell, we 
would have uJq = 1-37 meV, less than 2A2, and the ideal 
limit would amount to a good approximation. 



Sfqi is introduced as a scale factor for this situa- 
tion where an integral has been performed over q. The 
quantity lideaii'^) / S fqi is dimensionless. 



VIII. APPLICATION TO MAGNESIUM 
DIBORIDE 

For MgB2 the superconducting pairing comes primar- 
ily from optical phonons and is attractivcjii^ The intra 
a band effecting pairing potential D^^ is due to modes 
that become the doubly degenerate optical phonon of Eg 
symmetry at g w wherein the two B atoms in the unit 
cell move perpendicular to the c axis. The intra tt band 

(rt) 

pairing potential -D22 inter n — a band pairing po- 

tential are due to the mode that becomes the optical 
phonon of Big symmetry at g « wherein the B atoms 
move parallel to the c axis. Neither of these phonons 
for the relatively small (with respect to the inverse lat- 
tice constant) physical values of q are capable of mak- 
ing a contribution to the bubble correction, because the 
lines connecting the bubbles are, in fact proportional to 
phonon propagators and they must have the same sym- 
metry, Aig, as that of the bare Raman vertex. Thus £)'■'''' 
can only be due to the much weaker effects of Aig sym- 



metry acoustical phonons. We thus assume D 



IX. THE REAL CASE 

We now turn to the more realistic case, where Fj, Qj, 
and Pj have complicated dependences on lo and q. They 
depend on the properties of /(/3) = ^^^^'"^^ as a function 

Keeping in mind that 



of ^ = I3jiq,u,uj) 



4A5 



uj has an infinitesimally small positive imaginary part, 
we may write for various regions: 



in(yrT 



arcsin(/3) 



, when < 0, (45a) 



3in(yT^) 



fW) 



ln(/3 - 



when < /3 < 1, (45b) 

TTi/2 



1 



:, when /3 > 1. 

(45c) 



0. 



These expressions are to be inserted into Eqs. dM)) and 
(IT5t . We replace /3 with (2Aj)~ ^Jlo"^ — {qvju )2 and 
then average over u to obtain Fj{q,uj), Qj{q,uj), and 
Pj{q,uj). The u-integrals yield elliptic integrals for the 
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imaginary parts, but we know of no analytic forms for 
the real parts. We chose not to rely on numerical inte- 
gration to obtain the real parts of F, P, and Q. Instead 
we approximate the real parts of /, p, and r with func- 
tions of /3 that can be averaged over u. We used these 
the following approximate functions for the real part of 
/(/?). For /32 < 0, < /3 < ^, ^ < /3 < 1, and /3 > 1, 
respectively, we used: 



20 



fiO) (^) ^ 



7 



+ 



f^'^ (/?) 



3(1-/32) (i„^2)- 



P2 



1 

6^' 



(46a) 



(46b) 



(46c) 



(46d) 



These have the properties that /(°) {P) and /(^^ (/3) 
tend to 1 as /S^ tends to from below and above, re- 
spectively; /(2) (/3) has the same limiting, singular be- 
havior as /(/?), when /3 tends to 1 from below; and 
/(■^'' {13) tends to the correct limit, —1, when /3 tends 
to 1 from above. Moreover, the approximations /'■^•' (/3) 
and /(2) (^) agree when /? = 1/^/2. These approx- 
imations for /(/3) were inserted in Egs. p^ and 
to obtain approximate, integrable, versions of the real 
parts of p(/3) and r(/3). The averages over u were than 
calculated analytically with the help of 'Mathematica.' 
The resulting functions Fj , Pj , and Qj of q and w were 
used in Eas. (p4)) . ([27)) . and p5)) to calculate analytically 
Im[i3A(7).7* (<?, '^ + iO+)] —lra[x-f{q,oj + iQ'^)] that was in- 
serted into Eq. for numerical integration over q from 
to an upper limit of IQq" that was sufficiently large to 
be accurate in the peak region of the intensity plots. 
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FIG. 1: Scaled Raman intensity for the ideal case with inter- 
band pairing potential -D12 = 0.001, 0.003, and 0.01. The 
singularity at each peak was rounded. Also shown is the real 
case result for —D12 — 0.01. Here, a small positive imaginary 
part was added to the energy variable uj. 
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X. NUMERICAL RESULTS 

We start by comparing in Fig.([T]) the ideal case results 
for -D12 = 0.001, 0.003, and 0.01 with the real case 
result for the unrealistically small value of — D12 = 0.01. 

Note the similarity of the two —D12 — 0.01 spectra. 
Note also the discontinuity at w = 2A2 ~ 4.6 meV for 
the real case. For more realistic values of —D12 the effect 
of departure of the functions Fj , Pj , and Qj from their 
ideal case forms leads to dramatic changes in the Raman 
spectra, as shown in Fig. ([5]). Note that the discontinuity 
at a; = 2A2 is now joined by one at w = 2Ai = 13.5 
meV. 

Table I gives a calculated value for — Z3i2 of 0.08. 
The experimental result of of Blumberg et alM- shows 
that the peak occurs at 9.4 meV. Our calculation with 
—D12 = 0.08 gives the peak energy at 9.0 meV, whereas 
our calculation with — D12 = 0.10 has a peak energy very 
close to the experimental value. Taken together, the re- 
sults in Figs. (1-2) show how the collective mode, first 



FIG. 2: Scaled Raman intensity for the real case for —D12 = 
0.01, 0.02, 0.04, 0.08, 0.10, and 0.14. For the 0.01 plot a small 
positive imaginary part was added to the energy variable cj. 
'Wiggles' near 4.6 and 13.5 meV, the values of the 2 gaps, 
2Aiand 2A2, are artifacts of smoothing the numerical inte- 
gration carried out at a discreet number of points. 

predicted by Leggett, but corrected for final state inter- 
actions, evolves into a resonance mode peaked very close 
to the energy found by the Raman experiment as the pa- 
rameter —D12 increases from very small values to those 
close to the LDA result. 

We show a final comparison in Fig.(I5}D), namely be- 
tween the —D12 = 0.10 result from Fig. ([2]) and the result 
with the same value of —D12 obtained by setting q — 0. 
The effect of finite values of q on shifting up the energy 
of the peak is obvious. This is a kinematic effect that 
crudely results from the need for higher values of uj to 
compensate for nonzero values of {q ■ Vk)'^ in Eqs.([T]), 
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and p^ . 

These effects play out in an interesting way in our case 
where the 2 Fermi velocities differ by a factor of 12. The 
peak of -Im[x^(q, uj)] moves from 8 meV at g = to lower 
energy as q increases from to about 2 • 10'' m~^. Then 
the peak moves to higher energy. The behavior of the 
integrand in Eq. (jll[) can be seen in the density plot in 
Fig.dSi). 
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FIG. 3: (Color online) Results for the interband pairing po- 
tential —Di2 = 0.10. (a) Density plot of the integrand 
in Eq.(ll) with values of q (in units of 10^ m~^) extend- 
ing to lOg". Intensity scale is shown by the color bar on 
the right. The dashed curve starting at 4.6 meV obeys 
HvTrk = \JijS^ — The dashed slanted straight line obeys 

hv-nk = ijj, and the dashed curve starting at 13.5 meV obeys 
hvak = \/bj^ — 4A|. (b) The 'Integrated q' plot is the same 
as the — -Di2 = 0.10 plot in Fig. (2). The vertical dashed 
line at 9.5 meV in (a) denotes the value of Energy where the 
'Integrated q' plot has its maximum. The 'g = 0' plot is 
equivalent to a g = cut along the lower edge of the density 
plot in (a). 

Note that the shoulder at 12.3 meV in the 'Integrated 



q' curve of Fig ([SJd) corresponds to the peak at the q = 
lOg" cut at the top of the density plot of Fig. ([5^). If 
we were to extend the upper limit of the integral over q 
beyond lOg", this shoulder would move up in energy, but 
there would be no effect in the shape of the 'Integrated q' 
plot in the energy range between 4.6 and 11.5 meV. The 
dispersion of the peak in the density plot with wavevector 
q beyond about 10* m~^ is caused by the non-zero, but 
small, nature of the velocity Vc- If we had set = 0, the 
peak in the 'Integrated q' plot would occur at 9 meV, and 
it would be higher and narrower than in the plot shown 
in Fig.dnb). 

Recall that we had set D^^) = D^p) 7^ and D^'') = 
to obtain these results. Had we set D*^") =0 and D'^''^ = 
0, the resulting Raman intensity calculation would give 
a peak at a lower energy. In particular, the peak features 
shown in the the density plot in Fig.(|3^) would occur at 
energies 1.5-2 meV lower. Had we used D^^-' = D*^^^ = 
D^**) 7^ 0, they would have occured at still lower energies. 



XI. DISCUSSION 

To the author's knowledge this work amounts to the 
first attempt at a full treatment of finite wave vector ef- 
fects in the calculation of electronic Raman scattering 
in a superconductor. Our "real case" calculation takes 
these effects into account only through the use of the toy 
model assumption of spherical Fermi surfaces together 
with approximations to the real part of /(/3) [Eqs. (|46al 
I46dp ] . A more realistic calculation would involve a much 
more complex model of the Fermi surfaces and a great 
deal more massive numerical calculation. Nevertheless, 
the toy model used here seems to have captured the es- 
sential features of Leggett's mode as applied to the real 
multiband superconductor MgB2: damping due to decay 
into the n band(s) and finite wavevector effects in both 
TT and a bands. 

When applied to the normal state, our model does not 
account quantitatively or qualitatively for the experimen- 
tal results in Aig symmetryi^ The data show a rising con- 
tinuum in both Aig and E2g symmetries that is cut-off 
below 2A^ in the superconducting state. An explana- 
tion will require a model for the bare vertex 7^ that goes 
beyond the scope of the present work. 
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Appendix A: Flow, Counter-Flow, and the Screening 
Correction 

Here we derive an important result that has rather 
general applicability. It concerns the interplay between 
the correction for pairing and the screening correction 
when the wavevector q is zero. In that case we have 

. „x . „x . . „x arcsin(2f-) 




All the factors in square brackets in Eqs.Q, ([5]) and ([7]) 

become (w, 0) [T^ (t^, 0) - 2i/^k'i^-^Tf) (w, 0)] . 

To examine the implications for screening, we need to 
know the behavior of \ '{lu, 0), where the second sub- 
script denotes the value of 7^ = 1. We find that setting 

r\^l{w,0) = 2zAfcw"i and rj^^J(cj,0) = 1, 

thus making [T^,^\{u},0) - 2iAk'UJ-'^T^i^,\{uj,0)] = 0, triv- 
ially satisfies both Eq.(|4]) and [with the help of Eq.®] 
Eq. ([U . These results show that the term in square brack- 
ets vanishes for each k on the right side of Eq.© for 
Bj>xi{^,0)- Thus, i?^'ji(cj,0) = for any 7^. 

Of the two terms in square brackets in Eq.©, the T^, 
term can be thought of as 'flow' of a fc-dependent Ion- 
gitudinal particle/hole current, and the term as a 
fc-dependent particle/particle 'counter- flow' current pro- 
duced in response to the p/h current. For each value 
of k, p/p and p/h currents tend to cancel, and, when 
the p/h current is uniform in k space (7^. cx 1), we have 
just shown that this cancellation is exact for each k. At 
a deeper level the result F(^)(1) — 2iAuj~^ is a result of 
gauge invariance and particle conservation. For example, 
see the first of Nambu's Eqs. (7.6).— 

For non-zero q, Eqs. (H)), ^ and ([7]) apply with 
Pk{^,q) - fk{^,q) and rk{uj,q) - fk{i^,q) 7^ 0, and for 



small q, the leading terms in an expansion are oc g^. 
This has the consequence that the leading term in an 
expansion of i3^'_ri(w,q) is cx q^. This also applies to 
i3^*_ri(w, q), Bi^r-yi^,q), and _Bi.ri(i^, q) in the screen- 
ing correction. [See Eqs. (|9al and fT0)) .] Thus, the leading 
term in an expansion of the screening correction is cx . 

We have shown that the screening correction 
gives zero in the q — limit. We have sketched the 
proof here for the case of non-retarded interactions re- 
stricted to the Fermi surface. A similar proof involv- 
ing vertex corrections, bubble corrections, and screening 
should be possible for non-retarded interactions not re- 
stricted to the Fermi surface. Because of the underlying 
connection to gauge invariance, this statement about the 
screening correction should hold with retarded interac- 
tions as well. 

Thus, for Aig symmetry and g = 0, it is not cor- 
rect to use the bare bubble, Eq.©, B^i^^{uj,0) = 
—2 ^ dSk^'j^fk{^,Q)lki in the bare version of the screen- 
ing correction, Eq. (fTO| . with Fi — > I and F^ 7. One 
must have a solvable model for the pairing interaction 
[D^P\k,k')] and the vertex correction [£»(") (fc, fc')], and 
then solve equations like Eqs.© and ^ to obtain ex- 
pressions for F^^^(a;,0) and Fg(a;,0) to use in Eq.©. 
One way to do this would involve a re-interpretation 
of our solution for the multiband case. One can inter- 
pret matrix Eas. (|17ll?I|) as resulting from a fine-scale, 
or large n, discrete version of Eqs.(|4©. The solutions 
are given by Eqs (|22ll27p . The response function matri- 
ces R^P\uj,q), R'-P'''\u;,q), and R''P''''^\uj, q) are sym- 
metric. They can be shown to have the property that 
1 • R{uj, q) = R{u}, q) ■ 1 oc q"^ for small q. 

Another solvable model results from using the product 
ansatz D^^} cx AjAji. To obtain a solution, Eqs. ([19]) and 
(I20t must first be left-muhiplied by D^p'> = W^^-. 
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